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Abstract 

We calculate the excitation frequencies of the to = monopole and 
m = 2 quadrupole modes in the coUisionless regime by solving a non- 
linear Schrodinger equation for the condensate, coupled to a coUisionless 
Boltzmann equation for the quasiparticles. Since the dynamics of the non- 
condensate cloud is also taken into account, the theory satisfies the Kohn 
theorem. The spectrum turns out to be strongly temperature dependent 
and we compare our results with experiment. 

1 INTRODUCTION 

Since the experimental realization of a Bose condensate in trapped alkali va- 
pors, there has been renewed interest into the subject of degenerate quantum 
gases. Topics like the equilibrium properties of the condensate, the dynamics of 
condensate formation, topological defects, and collective excitations have been 
studied extensively. Of particular interest are the collective excitations in the 
coUisionless regime, because in this regime, where the mean free path of the Bo- 
goliubov quasiparticles is much larger than the wavelength of the collective exi- 
tations, there seems to be a discrepancy between the experimental observations 
and theoretical calculations. At temperatures far below the critical temperature 
Tc, measurements of the low-lying collective excitations ^ are in excellent 
agreement with theoretical calculations solving the Gross-Pitaevskii equation 
m, ^, m, which describes the condensate dynamics at zero temperature. At 
higher temperatures there is a considerable noncondensate fraction, and one 
has to include the mean-held interaction of the thermal cloud into the evolu- 
tion equation for the condensate wave function. Theoretical calculations solving 
the resulting nonlinear Schrodinger equation predict almost no temperature de- 
pendence of the lowest exitation frequencies ||^, ^, whereas experiments clearly 
show a large temperature dependence [^. This might partly be explained by 
including in the effective interaction between two colliding particles the many- 
body effect of the surrounding gas on the collisions, which causes the effective 
two-particle interaction to become strongly temperature dependent The 
frequencies of the low-lying modes will therefore also depend on temperature 
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pA[| . However, in these approaches the nonhnear Schrodinger equation describes 
the dynamics of the condensate in the presence of a static noncondensed cloud. 
As a result they violate the Kohn theorem, which states that there should al- 
ways be three center-of-mass modes with the trapping frequencies. Clearly, this 
violation is caused by the fact that we also have to describe the time evolution 
of the thermal cloud. Hence we propose to describe the collective excitations in 
the collisionless regime by a nonlinear Schrodinger equation for the condensate 
wave function that is coupled to a collisionless Boltzmann equation describing 
the dynamics of the noncondensed atoms |12, 13 1. This resolves our problem, 
because the resulting theory can be shown to contain the Kohn modes exactly. 
A full solution of the collisionless Boltzmann equation for the distribution func- 
tion of the Bogoliubov quasiparticles is rather complicated. Therefore we apply 
as a first step in this article the Hartree-Fock approximation for the quasipar- 
ticle dispersion. This is appropriate in the most interesting region, near the 
critical temperature Tc, where the mean-field interaction of the condensate is 
small compared to the average kinetic energy of the noncondensed cloud. We 
then determine the eigenfrequencies of the low-lying modes by a variational 
approach and compare these to the experiments. 



2 COLLISIONLESS DYNAMICS 

Our aim is to describe the coupled dynamics of the condensate and the ther- 
mal cloud in the collisionless regime, where the wavelength of the collective 
excitations is much larger than the mean free path of the quasiparticles. The 
time evolution of the condensate wave function <I>(x, t) describes the dynamics 
of both the condensate density no(x, t) = |$(x, and the superfluid velocity 
Vs(x,t) = n [$*(x,t)V$(x,t) - ^>(x,t)V^>*(x,i)] /2imno(x,t). In the mean- 
field approximation it obeys the nonlinear Schrodinger equation, that includes 
the effect of the mean-field interaction of the noncondensate density n'(x, t), 
and reads 



.^a^^ = |-^ + ^*™^(x)-^ + r2^[2n'(x,t)+no(x,t)]|cI.(x,a) 

Here, y*^'"P(x) denotes the external trapping potential. The factor of two differ- 
ence between the condensate and noncondensate mean-field contibutions results 
from the symmetrisation of the many-body wave function. For the nonconden- 
sate part this contributes a Hartree and a Fock term, but for the condensate 
part only a Hartree term. Moreover, the Hartree and the Fock contributions 
are equal because the two-body interaction is approximated by a hard-core 
potential y(x — x') = T^^5(x — x'), where the two-body scattering matrix 
rp2B _ 47r7i^a/m solves the Lipmann-Schwinger equation for the scattering of 
two particles with zero momentum. In the Thomas-Fermi limit one can neglect 
the average kinetic energy of the condensate relative to the mean-field inter- 
actions, and the equilibrium density profile of the condensate is approximately 
an inverted parabola. In the opposite limit where the mean-field interaction is 
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much less than the average kinetic energy, the result is essentially a gaussian 
density profile associated with the ground state of an harmonic oscillator. 

The Boltzmann equation is obtained by performing a gradient expansion on 
the equations of motion for the Wigner distribution, i.e. the Fourier transform 
of the one-particle density matrix. This is justified because the noncondensate 
density profile varies on a much larger length scale than the external trap- 
ping potential. It describes the time evolution of the quasiparticle distribution 
function F(k, x, f), which gives the noncondensate density profile, basically by 
integrating over momentum space. The local group velocity and the local force 
on a quasiparticle are given by the momentum and the spatial derivative of the 
dispersion, respectively. Therefore the Boltzmann equation reads 
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Sufficiently close to the critical temperature, where ksT ^ T^^no(x), we can 
treat the quasiparticles in the Hartree-Fock approximation and the dispersion 
is accurately given by 



u;(k,x,t) = ^ + y*™P(x)+2r2^[n'(x,t)+no(x,t)] . (3) 

The equilibrium distribution can be found by requiring the collision term to be 
equal to zero. For the equilibrium density profile the result is equivalent to a 
local density approximation. This is accurate for the conditions of interest, since 
also the coherence length of the gas is much smaller than the the scale on which 
the equilibrium profile varies. To find the collisionless exitation spectrum, we 
now consider fiuctuations around equilibrium. Since the eigenfrequencies of the 
lowest collisionless exitations are much larger than the average time between 
two collisions, we can neglect the contribution from the collisional term. The 
result is thus a linearization of the collionless Boltzmann equation, which is also 
known in the context of plasma physics as the Vlasov equation. 



3 VARIATIONAL APPROACH 



Variational calculations have been quite succesfull in obtaining the frequen- 
cies of the low-lying modes at zero temperature. Because of this, and because 
solving the coupled equations ([l|) and @j is rather difficult, we use a vari- 
ational method to solve these equations. Therefore we need two variational 
functions. One for the condensate wave function, the other for the quasiparti- 
cle distribution. Here, the external trapping potential is a quadratic potential 
ytrap^^^ _ '^.fjiij'^x'j /2. In the experiments performed at JILA and MIT the 
trap is cilindrically symmetric, hence uji = uj2 = oJr-, and oj^ = ujz. 

First, the condensate wave function is approximated by a gaussian, 



a>(x,i) 




hilSix'f 



(4) 
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This is appropriate for the experiments performed at JILA which are not in the 
Thomas-Fermi Hmit. In fact, the ratio of hu) and T^-^no is about one halve, 
where lj = (co^lUz)^^^ . However, choosing the condensate wave function to be a 
gaussian is known to give quantitatively good results for the condensate modes, 
even in the Thomas-Fermi limit where the ground state is well approximated 
by an inverted parabola. Therefore, we expect the above ansatz to be also at 
least qualitatively applicable to the MIT experiments. 

Second, to find an appropriate variational function for the quasiparticle 
distribution function, we take the distribution function to be an approximate 
Maxwell distribution for a gas of non-interacting bosons. To describe the 
monopole and quadrupole modes we introduce three scaling paramaters {ai} 
in the directions {xi} respectively, where i = {1,2,3}. Hence, 

F(k,x,i) = N' 

This results in a time dependent density profile n'{{xi/ai{t)}). For such a pro- 
file, the local current should be given by h{ki){'x.,t)/m = ai{t)xi/ai{t), which 
is correctly reproduced by our trial function. Furthermore, the equilibrium val- 
ues for Oi are in general different for i = {1.2,3} and larger than one, due 
to the interactions. This will cause the equilibruim density profile to broaden, 
but unfortunately also cause the distribution in momentum space to become 
anisotropic, which is incorrect because the mean-field interactions have no mo- 
mentum dependence. The constant overall factors {Cj} are introduced to com- 
pensate for this by choosing the value of Cj equal to 1/ a? g , where a^^e denotes 
the equilibrium value of ccj. 

We can now find evolution equations for the variational parameters {/?«, Aj, a^} 
by minimizing the lagrangian for the nonlinear Schrodinger equation and writ- 
ing the Boltzmann equation as a set of moment equations, that truncates be- 
cause of the variational ansatz. Our main result, the equations for {(3i},{Xi}, 
and {ai}, therefore reads 
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Here, Kth denotes the thermal wavelength {lirh/mkBT)^/'^ , lis equal to {lj.lz)^^^, 
li denotes the harmonic oscillator length {fi/muJiY/'^ , and Nq and A^' denote the 
total number of condensate and noncondensed atoms, respectively. 

To determine the spectrum of the low-lying collective excitations as a func- 
tion of r/Tc, we take the critical temperature and the number of condensate 
and noncondensed atoms to be given by the same expressions as those for an 
ideal gas in a trap, i.e. ksT^ = ^a)(iV/l. 202)^/3^ No = N [l - {T/T^f], and 
N' = N (T/TcY. In addition, we have included the effects of evaporative 
cooling, by approximating the dependence of the total number of particle N 
on T/Tc as found in the JILA experiment by a second order polynomial. We 
then find the fixed point of equations (6), (7) and (8) for each temperature and 
calculate the eigenvectors and their eigenvalues by linearizing around the fixed 
point. The results are plotted in figure 1. 
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Figure 1: The low-lying m = 0, 2 collisionless modes as a function of T/Tc, 
where <j) denotes the relative phase of density profiles of the condensed and 
noncondensed atoms. Also included are the experimental results for the m=0 
(triangles) and m=2 (circles) modes found in the JILA experiment P] 



4 DISCUSSION AND CONCLUSION 

Instead of the usual two modes, found when solving the Gross-Pitaevskii equa- 
tion or the nonlinear Schrodinger equation, there are now four modes. This 
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should not come as a surprise, since we are essentially dealing with two gas 
clouds. Without an interaction between these clouds, there would be a monopole 
m = and a quadrupole m = 2 mode for both the condensate and the noncon- 
densed cloud. If we turn on the interaction these modes get coupled, resulting 
in four modes where the condensate and the thermal cloud move either in or out 
of phase. These are the collisionless analogues of the hydrodynamic first and 
second sound modes [0, |l^ . Most importantly, we see from figure 1 that the 
temperature dependence of the m = 2 mode is in good agreement with the JILA 
experiment, in contrast with the previous approaches which do not take into 
account the dynamics of the noncondensed cloud. In addition, the m = mode 
does have the correct non-interacting limit near Tc, but the experimental data 
drops to the zero-temperature limit (10/3)^/^ ojr at a higher temperature than 
our theoretical curve, which nevertheless shows qualitatively the same behavior. 
We therefore believe that the avoided crossing between the in and out of phase 
monopole modes that causes this behavior, might be the reason for the strong 
temperature dependence of the m = mode found experimentally. The reason 
that here the quantitative agreement between our theory and the experiments 
is not that good, might be either the use of the Hartree-Fock approximation, 
or of the gaussian approximation for the condensate wave function. Work to 
improve on these approximations is in progress. Finally, we note that it might 
be possible to observe also the other two modes experimentally. Whether this is 
possible depends on the overlap of these modes with the applied perturbation, 
and on the damping of the modes, which we have neglected thusfar. However, 
in principle the collisionless Boltzmann equation also contains Landau damp- 
ing. Moreover, by including the collision term we should be able to describe 
collisional damping in the same way as Kavoulakis, Pethick and Smith jl6[ |. 
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